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Abstract: Quantum protocols often require the generation of specific quantum states. We 
describe a quantum algorithm for generating any prescribed quantum state. For an important 
subclass of states, including pure symmetric states, this algorithm is efficient. 

1 Introduction 

Many results in quantum information theory require the generation of specific quantum states, such as EPR 
pairs, or the implementation of specific quantum measurements, such as a von Neumann measurement in 
a Fourier transformed basis. Some states and measurements can be efficiently implemented using standard 
quantum computational primitives such as preparing a qubit in the state |0) and applying a sequence of 
quantum gates (from a finite set). EPR pairs can be prepared from the state |0)|0) using a Hadamard gate 
and a controlled-NOT gate. A von Neumann measurement in the Fourier basis can be efficiently realized by 
applying an inverse quantum Fourier transform and performing a von Neumann measurement in the standard 
computational basis, i.e. {|0}, |1)}. However many states and basis changes cannot be efficiently realized. This 
paper focusses on the generation of quantum states. For example, in [HMPEPC98], their improved frequency 
standard experiment requires the preparation of specific symmetric states on n qubits, where n is a parameter 
(number of ions). The algorithm we describe here will efficiently prepare the required symmetric state. This 
short paper will focus on the algorithm for generating the state, and will ignore issues related to errors 
and decoherence (for which the theory of fault-tolerant error-correction, or other stabilization methods, will 
apply). We do not have space to elaborate on the details of precision, but simple calculations that require 
0(log(i)) extra space and polylog(^) elementary operations allow us to generate any state with fidelity at 
least 1 — e. 

Suppose we want to generate the state |\&) = X^e{o 1}™ eVy " a x\x)- I n practice it suffices to generate a 
state | it) satisfying |($|$)| 2 > 1 — e for a given small real number e > 0. In this case, it suffices to 
approximate each ^y x and a x to accuracy poly(e) (i.e. 0(log(^)) bits of accuracy). Here we have factored 
out the phase in each term, and so the a x are all non-negative real values. Note that if we can prepare the 
state \Hf) = J2 X £{Q i}" a x\ x ), then we can approximate |\I/) arbitrarily well by introducing appropriate phase 
factors using methods discussed in [CEMM98]. We will therefore focus on a method for generating states 
|\I/) with non-negative real amplitudes. 

2 The algorithm 

In order to create the n-qubit state we will implement in sequence n controlled rotations, with the fcth 
rotation controlled by the state of the previous k — 1 qubits for k > 0. 

We will first define these controlled rotations, and then in the next section we will describe how we would 
implement them. 

Will extend the definition of a x to x £ {0, lp for 1 < j < n. Suppose we had a copy of and we 
measured the leftmost j qubits in the computational basis. Let a x be the non-negative real number so that 
a 2 equals the probability the measurement result is x. Then (a. XlX2 ... Xk _ 1 o/a XlX2 ... Xk _ 1 ) gives the conditional 
probability that the fcth qubit is |0), conditioned on the state of the first state of the k — 1 qubits being 
\x1X2 ■ ■ ■ Xk-i). Define a controlled rotation c — U XlX2 ^ Xk _ i0 by: 

\x 1 )\x 2 )...\x k ^ 1 )\o) c - u °>°>-°*->°., mm . . . 1^-0 ( "v?-?- 10 io) + a zz ii)) 
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As shown in Figure 1, the algorithm for generating the n-qubit state is a sequence of n such controlled 
rotations. It is easy to show by induction that after the first j controlled rotations are applied we have 
produced the state 
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\xiX 2 ■ ■ ■ Xj) 



XlK2---X 3 g{0,l}j 

and therefore after all n controlled rotations we have 



U* lX2 ... Xn _ l0 (c - t/JU...x n _ 2 o (• • • o- I7* (Itf |0» |0) . . .) |0)) |0) = |*). 
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Fig. 1. Network to generate 



3 Details 

In this section we show how to implement the controlled- f7* x with arbitrary precision. First assume we 

have a quantum register \^f) which encodes some "classical" description of the state VF The state |\&) must 
contain enough information to allow the probabilities o? x (or a related quantity, such as the cu x we define 
below) to be computed. We also use an ancilla register of 0(log(-)) qubits initialized to the state |0). Then 
we define operators Uk for each 1 < k < n as follows: 



where ojk satisfies 



|*)|0}|x!) . . . la*.!) ^ |*>|o; fc )|a;i> . . . \x k -i) 
cos 2 (27rw fc ) = (a xlX2 ... Xk _ l0 /a xlX2 ... Xk _ 1 ) 2 +0(poly(e)). 



A simple application of the techniques in [CEMM98] allow us to approximate (arbitrarily well) the transfor- 
mation: 

c - : |w)|Q) — > \u)e 2 ^ |0) , |w)|l) — > \u)e- 2 ^\l). 

" 



Also, define V 



. With these components, a network implementing c—U Xl Xk _ 1 o is shown in Figure 



2. Here we assumed that the algorithm works for a general family of states with classical descriptions |W). 
If we are only interested in producing a specific state |\&) the network can be simplified by removing the 
register containing |W) and simplifying each Uk to work only for that specific |^) (in the same way that one 
can simplify a circuit for adding variable inputs x and y to one that adds a fixed input 5 to variable input 

y)- 



4 Efficiency: an example 

The overall efficiency of our algorithm depends on how efficiently we can implement Uk', in other words, how 
efficiently we can compute the conditional probabilities (a XlX2 ... Xk o/a XlX2 ... Xk ) 2 or equivalently (a XlX2 ... Xk i/a XlX2 ... Xk Y 
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Fig. 2. Network implementing c— t/^ B 

One example for which this is easy is the symmetric states. The symmetric state \S r ) is defined to be an 
equally-weighted superposition of the computational basis states \x) that have Hamming weight H(x) = r 
(H(x) is the number of bits of x that equal 1). That is, 

\S r ) = —j= 

yj (™) H(x)=r 

The conditional probability ( flf2n^*^li ^ i s easily computed to be 

\ a »l*2- -*fc_l / 

r - H(xix 2 ■ ■ ■ Xk-i) 
n — k 

for 1 < k < n, and to be k — H(x\X2 ■ ■ ■ Xk-i) for k = n. The Hamming weight can be efficiently computed as 
shown in [KM01]. Then we simply need to reversibly compute the ui}. satisfying cos 2 (27ro;fc) = ( ^lf£n±^i£ ^ 
+ poly{e). 

Another example for which we can efficiently implement Uk is for more general symmetric pure states 

n 

where we are given the (3j values (as required in [HMPEPC98]). 

This technique will not allow us to generate efficiently all quantum states, but it will work for any family 
of states where for some reordering of the qubits we can efficiently compute the conditional probabilities 
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( a x 1 x 2 ...x k \ 
ttx 1 x 2 . -sc h J 
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